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1  Introduction 


It  is  well  established  that  thermal  energy  localization  resulting  from  the  rapid  deformation  of  het¬ 
erogeneous  reactive  solids,  such  as  propellants,  pyrotechnics,  and  high-explosives,  is  important  in 
the  ignition  and  onset  of  sustained  combustion  in  these  materials.  Mechanically  induced  ignition 
and  combustion  are  particularly  relevant  to  defense  related  applications  requiring  the  safe  handling 
and  storage  of  reactive  solids  and  the  development  of  modern  insensitive  munitions.  Localization 
typically  occurs  at  the  scale  of  heterogeneity  (i.e.,  meso-scale  or  grain  scale)  due  to  various  pro¬ 
cesses  including  inelastic  grain  and  binder  deformation,  intergranular  friction,  and  grain  fracture. 
Importantly,  energy  localization  at  the  meso-scale  can  have  a  significant  effect  on  the  large  scale 
(macro-scale)  material  response,  possibly  triggering  detonation  [16]. 

Mathematical  modeling  can  be  used  to  better  understand  the  interplay  of  localized  heating 
and  ignition  at  the  meso-scale  and  the  bulk  system  response  at  the  macro-scale.  Such  modeling 
necessarily  involves  the  coupling  of  physical  phenomena  occurring  over  disparate  length  and  time 
scales.  Several  bulk  models  have  been  developed  that  describe  mechanically  induced  transition 
to  detonation  for  both  pressed  and  granular  explosives  [1,  12,  15].  Often  these  models  are  only 
predictive  over  a  narrow  range  of  impact  conditions.  This  deficiency  is  largely  due  to  inaccurate 
descriptions  of  meso-scale  phenomena,  such  as  energy  localization  and  ignition,  also  referred  to  as 
hot-spot  formation,  within  the  context  of  the  bulk  models.  Some  models  give  a  more  complete 
description  than  others.  The  Baer  and  Nunziato  [1]  model  does  not  explicitly  account  for  hot-spot 
formation  but  specifies  an  ignition  criterion  that  is  based  solely  on  bulk  quantities.  The  model  of 
Johnson,  et  al  [12],  does  account  for  hot-spot  formation,  though  it  imposes  a  constant  value  for 
hot-spot  mass  fraction,  rather  than  describing  its  evolution,  and  estimates  hot-spot  temperature 
based  on  the  bulk  shock  pressure.  A  more  complete  description  is  given  by  the  model  of  Massoni, 
et  al  [15].  Their  composite  model  tracks  the  evolution  of  hot-spot  temperature  and  mass  fraction 
by  viscoplastic  pore  collapse,  in  a  manner  similar  to  that  first  proposed  by  Kang,  et  al  [14],  and 
couples  this  hot-spot  model  with  a  bulk,  two-phase  mixture  model.  However,  due  to  the  complexity 
of  their  composite  model,  it  is  unlikely  that  the  energetics  described  by  the  hot-spot  and  two-phase 
models  are  consistent.  Other  models  have  been  developed  that  track  growth  and  decay  of  hot-spots 
to  predict  explosion  times  and  thresholds  based  on  imposed  initial  conditions  for  hot-spot  size  and 
temperature  [3,  7,  23]. 

The  main  goal  of  this  study  is  to  develop  and  explore  a  technique  for  formulating  a  comprehen¬ 
sive  and  predictive  model  of  mechanically  induced  hot-spot  formation  and  ignition  in  heterogeneous 
reactive  solids.  The  technique  requires  1)  a  model  for  the  bulk  material  response,  2)  a  model  for  the 
meso-scale  structure  (e.g.,  grain  size,  grain  packing,  etc.),  3)  a  localization  strategy  for  depositing 
bulk  dissipated  mechanical  energy  at  the  meso-scale,  and  4)  a  model  for  the  meso-scale  response. 
The  meso-scale  response  will  allow  for  the  evolution  of  hot-spot  temperature  and  mass  fraction, 
subject  to  the  localization  strategy,  rather  than  imposing  their  values.  Here,  it  is  assumed  that 
the  bulk  material  response  is  experimentally  well-characterized  and  is  accurately  predicted  by  the 
bulk  model;  as  such,  we  will  maintain  the  integrity  of  bulk  model  predictions  and  require  that  the 
integrated  mass,  momentum,  and  energy  of  the  meso-scale  model  locally  equals  that  given  by  the 
bulk  model.  Moreover,  this  constraint  is  consistent  with  the  common  interpretation  that  the  bulk 
response  is  an  average  manifestation  of  the  meso-scale  response.  Models  that  do  not  explicitly 
enforce  this  constraint  may  give  inconsistent  meso-scale  and  bulk  predictions.  While  such  a  con¬ 
straint  may  be  overly  restrictive,  it  does  provide  a  rational  framework  from  which  to  formulate  and 
explore  the  implications  of  various  localization  strategies. 
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In  this  study,  some  rigor  is  sacrificed  in  favor  of  tractability.  Rather  than  considering  a  heteroge¬ 
neous  system  of  explosive  grains  embedded  within  a  plastic-like  binder  (e.g.,  PBXs),  we  focus  only 
on  the  explosive  grains.  The  response  of  granular  explosives  to  dynamic  loading  is  of  interests  as 
they  are  typically  the  main  load  bearing  constituent  of  PBXs.  In  particular,  we  model  the  dynamic 
compaction  of  granular  HMX  (cyclo-tetramethylene-tetranitramine)  because  it  is  both  commonly 
used  in  aerospace  and  defense  related  applications  and  its  quaisi-static  compaction  behavior  is  exper¬ 
imentally  well-characterized  [5].  Further,  the  dynamic  compaction  and  Deflagration-to- Detonation 
Tra.Tisit.inn  (DDT)  of  granular  HMX  have  been  extensively  studied  [1,  2,  4,  9,  16,  18,  21].  Experi¬ 
mental  data  indicate  that  sustained  combustion  and  DDT  of  confined  granular  HMX  can  be  readily 
triggered  by  compaction  waves  resulting  from  weak  mechanical  impact  (piston  speeds  ~  80—90  m/s) 
[16,  22];  thus,  a  fairly  well-defined  threshold  exists  for  comparing  our  model  predictions  against. 
Other  modeling  simplifications  made  in  this  study  are  discussed  in  appropriate  sections  of  this 
report.  Though  this  study  focuses  only  on  granular  HMX,  the  technique  can  be  generalized  and 
applied  to  systems  containing  a  binder  phase  as  the  physical,  modeling,  and  computational  issues 
involved  are  similar. 

An  outline  of  this  report  is  as  follows.  First,  a  brief  discussion  of  dynamic  compaction  is 
given  and  key  issues  related  to  thermal  energy  localization  and  ignition  are  identified.  Second, 
the  mathematical  model  is  posed  and  the  localization  strategy  explained.  Third,  the  model  is 
numerically  solved  to  give  predictions  for  localized  heating  and  ignition  of  granular  HMX.  The 
piston  speed  needed  to  induce  sustained  combustion  is  predicted  and  compared  to  experiments, 
the  sensitivity  of  the  model  to  variations  in  key  localization  parameters  is  demonstrated,  and  the 
influence  of  grain  size  on  the  predicted  threshold  for  sustained  combustion  is  determined.  Lastly, 
some  conclusions  and  suggestions  for  future  work  are  given. 

2  Dynamic  Compaction 

A  common  paradigm  for  the  dynamic  compaction  of  a  granular  material  is  the  piston  impact 
problem  shown  in  Fig.  1(a).  In  the  figure,  the  uncompacted  material  has  30%  porosity  {4'o  =  0.70) 
and  is  in  an  unstressed  state;  this  initial  porosity  is  representative  of  loose  granular  energetic 
solids.  Compaction  is  induced  by  a  piston  moving  with  constant  speed  {Up  =  100  m/s).  Upon 
piston  impact,  compaction  energy  is  transmitted  through  the  material  by  shear  and  compression 
waves  that  involve  interactions  between  grains  through  their  contact  surfaces,  as  illustrated  in 
Fig.  1(b).  These  local  grain  interactions  result  in  both  stress  and  temperature  fluctuations  within 
the  material  meso-scale  structure.  The  integrated  effect  of  these  fluctuations  gives  rise  to  a  bulk 
wave  that  propagates  away  from  the  piston  {D  «  400  m/s)  leaving  the  material  in  a  stressed, 
compacted  state  with  low  porosity  {4>c  «  0.96);  this  wave  is  commonly  referred  to  as  a  dynamic 
compaction  wave.  The  compaction  process  is  not  instantaneous,  but  occurs  over  a  finite  length 
zone  (d  4  mm)  referred  to  as  the  compaction  zone.  Within  the  compaction  zone,  the  granular 
material  is  rapidly  deformed  and  heated  due  to  dissipative  interactions  between  grains;  common 
heating  mechanisms  include  intergranular  friction,  grain  fracture,  plastic  deformation  of  grains,  and 
compression  of  grains.  Temperature  fluctuations  within  the  compaction  zone  can  be  significant.  If 
the  deformation  rates  are  sufficiently  high,  these  dissipative  processes  will  lead  to  the  formation  of 
thin,  localized  regions  of  intense  heating  near  the  contact  surfaces  between  grains  as  there  is  little 
time  for  thermal  cooling  by  conduction  into  the  cooler  grain  interiors.  For  example,  using  the  values 
depicted  in  Fig.  1,  the  bulk  strain  rate  of  the  granular  material  induced  by  the  compaction  wave  is 
approximately  given  by  e  =  (1  -  ^o/</’c)/T  =  0.203  where  T  =  d/{D  -  Up)  =  13.3  )us  is  the 
residence  time  of  a  grain  within  the  compaction  zone.  It  will  be  shown  in  Section  4  that  such  strain 
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Figure  1:  Schematic  of  a  dynamic  compaction  process:  (a)  propagation  of  a  bulk  compaction  wave; 
(b)  localized  stresses  and  grain  heating  at  the  contact  surface  between  grains. 


rates  can  result  in  significant  volumetric  heating  rates  in  the  vicinity  of  intergranular  contact  points 
(e.g,,  5.0  X  10^  GW/m^).  Using  the  thermal  diffusivity  of  solid  HMX,  a  =  1.76  x  10“^  m^/s,  and 
assuming  that  all  heating  occurs  at  intergranular  contact  surfaces,  the  distance  that  thermal  energy 
is  conducted  into  a  grain  within  the  compaction  zone,  identified  by  S  in  Fig.  1(b),  is  approximately 
S  =  -s/aJT  ===  1.5  /xm  <  i,  where  I  ^  100  yum  is  an  average  grain  size.  Consequently,  for  the  case 
shown  here,  thermal  energy  will  be  localized  near  grain  contact  surfaces  within  the  compaction  zone. 
The  existence  of  these  fine-scale  features  poses  significant  modeling  and  computational  challenges 
that  will  be  described  later  in  this  report.  Compaction  induced  thermal  energy  localization  in 
energetic  solids  can  be  a  major  mechanism  for  combustion  initiation.  If  a  sufficiently  large  number 
of  intense  localization  sites  exist,  prompt  initiation  of  sustained  combustion  will  result  and  explosion 
is  possible. 

3  Mathematical  Model 

A  mathematical  model  for  compaction  induced  hot-spot  formation  is  outlined  in  this  section.  The 
bulk  compaction  model  is  first  given  followed  by  the  meso-scale  energy  localization  model.  Through¬ 
out  this  section,  specific  functional  dependencies  of  some  variables  are  emphasized  to  highlight  the 
coupling  of  the  bulk  and  meso-scale  models. 

3,1  Bulk  Compaction  Model 

A  one-dimensional  (1-D)  bulk  compaction  model  is  used  in  this  study.  Though  there  can  exist 
substantial  stress  bridging  between  grains  at  the  meso-scale  [17,  20],  resulting  in  spatially  3-D  flow, 
the  spatially  averaged  response  can  be  approximated  as  1-D  provided  that  appropriate  boundary 
conditions  are  specified.  Stress  bridging  is,  however,  an  important  issue  concerning  energy  local¬ 
ization  that  will  be  discussed  later.  The  bulk  compaction  model  was  formulated  by  Gonthier  et  al 
[9];  as  such,  the  reader  is  referred  there  for  details.  Here,  we  briefly  summarize  the  bulk  model, 
focusing  on  issues  relevant  to  compaction  induced  localized  heating. 
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The  bulk  model  equations  are  given  by 


d 

(  p  ^ 

dt 

pu 

[pE  J 

dx 

\ 

fm 

pu^-\-P 


\ 


=  0, 


d(j>  d(l) 
dt 


lie 


0 


if  /  > 

otherwise. 


(1) 

(2) 

(3) 


Independent  variables  these  equations  are  time  t  and  position  x.  Dependent  variables  for  the 
granular  solid  are  density  p;  velocity  u;  pressure  P;  total  mass  specific  energy  E  =  e  +  v? 
where  e  is  the  mass  specific  internal  energy;  solid  volume  fraction  no-load  volume  fraction 
intragranular  stress  and  equilibrium  no-load  volume  fraction  /(<?!>).  Quantities  associated 

with  the  pure  phase  solid,  denoted  by  subscript  “s”,  are  related  to  the  granular  solid  variables 
by  ps  ^  p/(j),  Ps  ^  P/^,  and  Cs  e-  B,  where  P(^  -  <^)  =  f#'  is  recoverable  compaction 

energy.  The  parameters  pc  and  p  appearing  in  Eqs.  (2)  and  (3)  govern  the  rates  of  relaxation  to  the 
equilibria  Ps  =  ^  and  /  =  respectively.  The  system  of  equations  is  mathematically  closed  given 
equations  of  state  Ps(ps,T)  and  es(ps,T),  where  T  is  the  bulk  temperature  of  the  granular  solid. 
In  this  study,  we  slightly  modify  the  inert  compaction  model  of  Ref.  [9]  to  account  for  chemical 
reaction  and  take  e5(T,  A)  =  c^T  -  Ag,  assuming  an  incompressible  solid,  where  Cy  is  the  specific 
heat  at  constant  volume,  A  is  a  bulk  reaction  progress  variable  that  is  determined  by  the  meso- 
scale  response,  and  q  is  the  specific  heat  of  reaction;  the  inclusion  of  chemical  reaction  is  further 
explained  in  the  following  section. 


Equation  (1)  is  a  system  of  conservation  expressions  for  the  granular  solid  mass,  momentum, 
and  total  energy.  Equations  (2)  and  (3)  are  evolution  equations  for  (j)  and  constructed  so  that 
irreversible  compaction  is  compatible  with  the  Second  Law  of  Thermodynamics.  It  is  noted  that  for 
an  incompressible  solid,  the  mechanical  response  of  the  granular  material,  given  by  the  mass  and 
momentum  expressions  of  Eq.  (1),  and  Eqs.  (2)  and  (3),  is  independent  of  the  thermal  response. 
This  note  will  prove  useful  later. 


We  focus  on  compaction  energetics  which  is  important  for  the  development  of  the  energy  local¬ 
ization  model.  Because  low  pressure  compaction  (<  100  MPa)  associated  with  weak  initiation  of 
DDT  induces  only  small  density  changes,  we  assume  an  incompressible  solid  and  ignore  compres¬ 
sion  work.  In  this  case,  it  can  be  shown  from  Eqs.  (1-3)  that  the  evolution  of  internal  energy  for 
the  granular  material  is  given  by 

—  —  ^  j_  —  (4) 

dt  dt  dt  ’ 

where 

dcs  ^  jPs  "  /?)  #  ^d(l>  .g. 

dt  p  dt  p  dt' 

(6) 

dt  pdt 

Here,  d/dt  —  dldt  +  ud/dx  is  the  convective  derivative.  Equations  (5)  and  (6)  describe  the  evolu¬ 
tion  of  bulk  thermal  energy  and  recoverable  compaction  energy,  respectively.  We  are  particularly 
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(c)  (d) 

Figure  2:  Predicted  variation  in  bulk  properties  through  the  compaction  zone  for  the  inert  com¬ 
paction  of  granular  HMX  by  a  Up  =  85.6  m/s  piston:  (a)  solid  volume  fraction,  (b)  velocity,  (c) 
pressure,  and  (d)  temperature. 

interested  in  Eq.  (5)  as  bulk  thermal  energy  will  be  deposited  at  the  meso-scale  based  on  our  lo¬ 
calization  strategy.  Entropy  considerations  indicate  that  each  term  on  the  right  side  of  Eq.  (5)  is 
dissipative. 

To  illustrate  features  of  the  bulk  compaction  model,  representative  predictions  for  the  spatial 
variation  in  the  bulk  quantities  u,  Pg,  and  T  through  the  compaction  zone  of  an  inert,  steady, 
dynamic  compaction  wave  in  granular  HMX  (A  =  0)  are  given  in  Fig.  2  for  the  piston  speed 
Up  =  85.6  m/s.  The  numerical  solution  technique  used  to  solve  the  model  equations  is  discussed 
in  Section  4  of  this  report.  Here,  the  compaction  wave  is  shown  propagating  to  the  right  at  speed 
D  =  297.5  m/s.  A  continuous,  dispersed  wave  structure  is  predicted  due  to  compaction  induced 
mechanical  dissipation.  Similar  low  speed,  dispersed  structures  have  been  experimentally  observed 
[16,  21]  and  predicted  by  detailed  meso-scale  numerical  simulations  [17]  for  granular  HMX.  Smooth 
increases  from  their  initial  values  are  predicted  for  all  quantities  through  the  compaction  zone.  The 
solid  volume  fraction  increases  from  its  initial  free  pour  value  of  =  0.655  to  a  final  equilibrium 
compacted  value  of  (/>  =  0.92.  The  velocity  increases  from  zero  to  Up  as  required  by  the  rear 
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boundary  condition  at  the  piston  surface.  The  initially  unstressed  granular  material  reaches  an 
equilibrium  stress  of  Ps  =  34.5  MPa  in  the  compacted  state.  Importantly,  the  bulk  temperature 
rise  associated  with  the  compaction  process  is  only  about  AT  =  2.4  K  which  is  much  too  low 
to  initiate  chemical  reaction.  Further,  the  bulk  temperature  rise  for  a  piston  speed  of  400  m/s 
is  only  approximately  50  K.  As  such,  combustion  models  that  are  based  on  bulk  temperature 
are  too  insensitive  to  accurately  predict  the  onset  of  self-sustained  combustion  for  varying  impact 
conditions.  This  issue  motivates  the  need  to  develop  rational  energy  localization  models.  It  is 
shown  in  Ref.  [9]  that  predictions  of  the  bulk  compaction  model  for  equilibrium  end  states  and 
compaction  wave  speeds  agree  well  with  experimental  data  for  a  wide  range  of  piston  speeds. 

3.2  Energy  Localization  Model 

In  this  subsection  we  describe  a  thermal  energy  localization  model  that  is  consistent  with  the  bulk 
model.  The  localization  model  consists  of  1)  a  model  for  the  evolution  of  meso-scale  structure,  2)  a 
strategy  for  localizing  bulk  dissipated  energy  at  the  meso-scale,  and  3)  a  model  for  predicting  the 
meso-scale  response. 


3.2.1  Meso~Scale  Structure 


The  meso-scale  structure  of  granular  HMX,  and  PBXs  containing  HMX  and  binder,  consists  of  a 
broad  distribution  of  grain  sizes  and  shapes  packed  in  complex  arrangements.  Experiments  indicate 
that  substantial  changes  in  the  meso-scale  structure  occur  as  the  material  is  mechanically  loaded 
[4,  19].  These  changes  are  in  large  part  due  to  the  brittle  nature  of  HMX  that  readily  fractures 
when  compressed;  grain/binder  deconsolidation  is  also  observed  in  PBXs.  A  detailed  model  for  the 
evolution  of  grain  size  distribution  that  accounts  for  grain  fracture  can  be  formulated  using  integral 
techniques  [8],  but  the  models  are  complex,  computationally  expensive  to  solve,  and  beyond  the 
scope  of  this  study.  For  now,  we  ignore  grain  fracture  and  track  the  evolution  of  grain  number 
density  n{x^t)  by 

|  +  ^(n.,=0.  (7) 


Fracture  can  be  accounted  for,  within  the  context  of  a  size  distribution  function,  by  supplying  a 
forcing  term  F dm  in  Eq.  (7),  where  F{m,x,t)  is  the  volumetric  fracture  rate  per  unit  grain 
mass  that  must  satisfy  mPdm  =  0  for  mass  conservation;  this  topic  will  be  the  focus  of  future 
work.  We  additionally  assume  that  the  grains  are  incompressible  and  of  uniform  size  and  spherical 
shape;  they  can  be  characterized  by  their  constant  radii  R.  Grain  number  density  is  related  to 
solid  volume  fraction  and  grain  size  by  n  =  ^/(4/37ri?^).  We  do  not  give  a  detailed  description 
of  the  grain  packing  arrangement,  but  specify  the  average  number  of  intergranular  contact  points 
(surfaces)  per  grain,  7.  As  an  illustration,  a  face  centered  cubic  (FCC)  arrangement  of  uniform 
size  spherical  grains  has  7  =  12,  whereas  a  simple  cubic  (SC)  arrangement  has  7  =  6.  Though 
a  real  granular  material  will  not  have  such  an  ordered  arrangement,  it  is  reasonable  to  assume 
1  <  7  <  16.  In  summary,  the  meso-scale  structure  is  characterized  in  this  study  by  the  variable  (j) 
and  the  constant  parameters  R  and  7. 


3.2.2  Localization  Strategy 

An  appropriate  energy  localization  strategy  should  include  information  based  on  the  dissipative 
mechanics  occurring  at  the  material  meso-scale.  Applied  bulk  loads  are  transmitted  at  the  meso- 
scale  by  intergranular  contact.  Experiments  [20]  and  numerical  simulations  [17]  indicate  that  the 
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Property 

Value 

£■,  Young’s  Modulus 

24.0  GPa 

y,  Yield  Strength 

0.37  GPa 

Poisson’s  Ratio 

0.2 

ps,  Density 

1903.0  kg/m^ 

Cy,  Specific  Heat 

1.5  kJ/kg/K 

A:,  Thermal  Conductivity 

0.502  W/m/K 

Table  1:  Material  properties  for  HMX. 


bulk  load  is  not  uniformly  transmitted  through  the  meso-structure,  but  is  preferentially  transmit¬ 
ted  along  certain  paths  that  form  stress  chains;  grains  not  involved  in  these  chains  remain  mostly 
unstressed  and  inactive.  It  is  reasonable  then  to  expect  that  grain  size,  shape,  and  packing  arrange¬ 
ment  may  have  a  large  effect  on  the  meso-scale  stress  state.  The  stress  state  will  also  depend  on  the 
magnitude  of  the  applied  load,  and  possibly  the  loading  rate.  For  weak  loads,  the  meso-structural 
response  will  mostly  involve  localized  elastic  deformation  and  friction  near  active  intergranular 
contact  points  (surfaces).  In  this  case,  the  elastic  stress  field  in  the  vicinity  of  the  contact  surface 
is  approximated  by  classical  Hertz  contact  theory  [13].  The  onset  of  plastic  deformation  will  occur 
within  the  grain,  near  the  contact  surface,  at  the  location  where  the  principal  shear  stress,  r ,  first 
equals  the  material  yield  strength,  Y.  An  expression  for  t  based  on  the  Hertz  solution  is  given  by 


T  =  -Pc 
2 

(1  +  ^^) 

[’-( 

9  (t). 

(8) 


where 


a  = 


nR*Pc 
2E*  ’ 


(9) 


and 


E*  = 


E 

2(1-  1/2)  ■ 


Here,  E  is  Young’s  modulus,  v  is  Poisson’s  ratio.  Pc  is  the  stress  at  the  contact  center,  and  a; 
is  distance  into  the  grain  normal  to  the  planar,  circular  contact  surface  of  radius  a  as  shown  in 
Fig.  3(a).  We  have  assumed  that  contact  occurs  between  spheres  of  equal  size.  For  axi-symmetric 
contact,  the  onset  of  plastic  deformation  occurs  when  Pc  =  1.6y,  which  is  based  on  von  Mises’  yield 
criterion.  Both  the  maximum  shear  stress  and  its  location  within  the  grain  at  the  onset  of  plastic 
deformation  can  be  determined  from  Eq.  (8).  Using  the  material  properties  for  HMX  listed  in  Table 
1,  the  onset  of  plastic  deformation  for  a  grain  size  of  i?  =  25  fj,m  will  occur,  for  Pc  =  0.59  GPa  at  a 
depth  z  =  0.45a  =  0.93  jum  corresponding  to  (r)max  =  0.33Pc  =  0.2  GPa.  Continued  loading  will 
cause  the  volume  of  plastically  deformed  material  to  increase  resulting  in  plastic  heating.  Grain 
fracture  and  intergranular  friction  may  also  be  significant,  particularly  at  elevated  loads.  Bulk 
dissipated  mechanical  energy  is  the  integrated  effect  of  these  meso-scale  dissipative  processes. 


A  phenomenological  localization  model  can  be  developed  based  on  the  foregoing  discussion.  In 
so  doing,  we  will  make  several  simplifications  and  explore  their  consequences.  Because  both  plastic 
deformation  and  friction  occur  in  the  neighborhood  of  intergranular  contact  points  (surfaces),  we 
will  track  solid  regions  surrounding  these  points  within  which  bulk  dissipated  mechanical  energy 
will  be  thermalized;  these  regions  will  be  referred  to  as  energy  localization  centers.  The  number  of 
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localization  centers  per  unit  volume,  ndx.t),  is  related  to  the  number  of  contact  points  per  grain 
and  grain  number  density  by  =  ^7^?  the  prefactor  1/2  is  introduced  because  each  localization 
center  involves  contact  between  two  grains.  The  localization  centers  are  evenly  distributed  and 
individually  surrounded  by  solid  spheres  of  radii  ro,  as  shown  in  Fig.  3,  where  Vq  =  The 

above  expression  for  Uc  can  be  combined  with  this  expression  for  Tq  to  obtain  ric  =  thus, 

all  solid  mass  is  encompassed  by  the  localization  spheres.  For  an  incompressible  solid,  as  assumed  in 
this  study  (i.e.,  R  =  constant),  is  constant  for  a  given  value  of  7.  Because  grain  contact  surfaces 
in  real  materials  would  not  be  evenly  distributed  as  assumed  here  due  to  stress  chain  formation 
and  grain  size  variations,  bulk  thermal  energy  would  be  localized  at  a  fewer  number  of  closely 
packed  localization  sites;  there  may  also  exist  significant  thermal  interaction  between  these  sites. 
The  present  model  can  indirectly  account  for  increased  localization  by  reducing  the  value  of  7,  but 
does  not  account  for  asymmetrical  thermal  interaction  between  sites  as  discussed  in  the  following 
subsection.  We  note,  however,  that  a  comparison  of  ignition  predictions  with  experimental  data  for 
HMX  suggests  that  the  onset  of  sustained  bulk  combustion  may  be  reasonably  insensitive  to  the 
spatial  distribution  of  localization  sites.  Nonetheless,  the  model  can  be  generalized  in  the  future  to 
more  accurately  reflect  these  features. 


Thermal  energy  will  be  deposited  over  a  volume  of  radius  ^  ^0  which  defines  the  local¬ 

ization  center.  The  initial  value  for  Vc  is  taken  as  the  radius  of  the  intergranular  contact  surface,  a, 
at  the  onset  of  plastic  deformation  within  the  grain.  This  volume  is  close  to  the  elliptical  volume 
defined  by  a  and  the  distance  from  the  contact  surface  at  which  plastic  deformation  first  occurs, 
=  0.45a.  Thus,  we  have  from  Eq.  (9)  that 


rdx,0) 


ttR^Pc 
2E*  ’ 


(10) 


where  Pc  —  1.6y.  This  assumption  is  reasonable  in  that  prior  to  the  onset  of  plastic  deformation 
most  dissipated  energy  will  be  due  to  intergranular  friction  and  will,  thus,  be  localized  near  the  con¬ 
tact  surface  within  the  region  r  <  rdx,0).  It  is  noted  that  the  initial  localization  volume  decreases 
with  grain  size.  To  obtain  an  expression  for  the  evolution  of  Vc  due  to  material  compaction,  we 
assume  that  the  material  is  perfectly  plastic  and  is  in  an  uncontained  mode  of  plastic  deformation; 
as  such,  details  of  the  transitional  range  of  loading  associated  with  purely  elastic,  elastic-plastic, 
and  fully  plastic  flow  are  not  modeled.  The  plastic  flow  stress  is  taken  to  be  Py  =  3.0y  [13].  This 
expression  for  flow  stress  has  also  been  used  to  describe  the  quasi-static  compaction  behavior  of 


Figure  3:  Schematic  of  the  thermal  energy  localization  strategy. 
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ball  powder  [11].  We  equate  the  volumetric  rate  of  work  done  by  the  plastic  flow  stress  to  the  bulk 
volumetric  dissipated  energy  given  by  Eq.  (5): 


n  I 

’'cPv-jf  = 


(11) 


where  Vc  =  |7rr^  is  the  localization  (plastic)  volume  for  a  single  sphere.  An  expression  for  the 
evolution  of  Vc  is  thus  given  by 

dvc _  Ps  ^  dcs  fl2) 

dt  A'Kr^^UcPY  dt  ’ 

where  Py  =  3.0y.  Though  we  have  assumed  that  all  bulk  dissipated  energy  is  the  result  of 
plastic  deformation,  the  bulk  dissipated  energy  can  be  easily  partitioned  into  both  frictional  and 
plastic  components;  this  would  involve  minimal  modification  of  the  model.  The  model  sensitivity 
to  variations  in  both  the  initial  localization  volume  and  the  plastic  flow  stress  is  addressed  in 
Subsection  4.3. 


3,2.3  Meso-Scale  Response 

A  model  for  the  meso-scale  response  to  localized  heating  will  now  be  described.  The  approach  is 
similar  to  the  modified  Prank-Kamenetskii  thermal  explosion  theory  of  Foster,  et  al  [6].  To  this  end, 
we  track  the  time-dependent,  radial  evolution  of  thermal  energy  and  chemical  reaction  progress  of 
HMX  contained  within  a  single  localization  sphere  subject  to  thermal  energy  deposition  based  on 
the  strategy  of  Subsection  3.2.2.  Evolution  equations  for  thermal  energy  and  reaction  progress  are 
given  by  the  following,  respectively: 


Here,  r  is  radial  position  within  the  localization  sphere  and  djdt  =  djdt  -Vudjdx  is  the  convective 
derivative.  It  is  necessary  to  include  the  convective  derivative  because  the  grains,  and  thus  the 
localization  spheres,  locally  propagate  at  the  bulk  velocity  u.  Variables  labeled  with  a  “hat”  (•) 
are  associated  with  the  localization  sphere  and  vary  not  only  with  x  and  t  but  also  with  r;  they 
include  solid  temperature,  T{x^  r,  t),  reaction  progress,  A(a:,  r,  t)  (0  <  A  <  1,  where  A  =  1  is  complete 
reaction),  and  volumetric  thermal  energy  deposition  rate  S{x^  r,  t).  The  meso-scale  and  bulk  models 
will  be  coupled  through  the  energy  source  term  S.  Constant  parameters  contained  in  Eqs.  (13) 
and  (14),  and  not  yet  defined,  are  thermal  conductivity,  Arrhenius  prefactor,  z,  and  activation 
temperature,  Ta.  Equation  (13)  is  consistent  with  the  caloric  equation  of  state  =  CyT  -  Xq, 
Equation  (14)  models  a  one-step,  irreversible  chemical  reaction  having  Arrhenius  kinetics.  Values 
for  the  thermal  and  chemical  parameters  for  HMX  are  given  in  Tables  1  and  2,  respectively. 

To  maintain  consistency  between  the  meso-scale  and  bulk  models,  we  require  that  the  evolution 
of  mass,  linear  momentum,  and  thermal  energy  at  the  meso-scale  locally  equals  that  given  by  the 
bulk  model.  It  can  be  easily  shown  that  the  mass  constraint  is  identically  satisfied  based  on  the 
expression  for  Vq  given  in  Subsection  3.2.2.  Also,  because  we  are  primarily  interested  in  the  meso- 
scale  thermal  energy  response,  we  have  assumed  that  all  grains  move  at  the  local  bulk  velocity; 
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Property 

Value 

Arrhenius  Prefactor 

5.00  X  lOi^  s-i 

q,  Specific  Heat  of  Reaction 

5.95  X  10*^  J/kg 

Tfl,  Activation  Temperature 

2.65  X  10^  K 

Table  2:  Chemical  properties  for  HMX. 


consequently,  the  linear  momentum  constraint  is  trivially  satisfied.  It  remains  to  satisfy  the  thermal 
energy  constraint  given  by 

(^J  Ps'^c  47r  J  drdc^  .  (15) 

The  left  hand  side  of  this  equation  is  the  evolution  of  bulk  thermal  energy  for  a  volume  element  of 
arbitrary  length  in  the  a:-direction,  where  Cg  =  CyT  —  \q  is  the  bulk  specific  internal  energy.  The 
right  hand  side  of  this  equation  is  the  evolution  of  integrated  thermal  energy  at  the  meso-scale. 
Using  the  relation  Uc  =  Eq.  (15)  gives  by  inspection 

3 

es{x,t)  =  -o  r^es{x,r,t)  dr;  (16) 

r$  Jo 

thus,  for  an  incompressible  solid,  the  bulk  internal  energy  is  the  volumetric  average  of  the  energy 
contained  within  a  localization  sphere.  Likewise,  the  following  relations  for  the  bulk  temperature 
and  reaction  progress  variable  can  be  obtained  from  Eq.  (16)  using  the  expressions  for  and  eg'- 


3  /■’■<’ 

T(x,t)  =  4/ 
rt  Jo 

3  ro 

=  4 


r^f{x^r,t)  dr, 


r^A(a?,r,t)  dr. 


(17) 

(18) 


The  expression  for  S{x,  r,  t)  must  be  chosen  so  that  Eq.  (16)  is  identically  satisfied.  To  this  end, 
we  first  take  the  convective  derivative  of  Eq.  (16),  with  constant,  to  obtain 


dcg 

dt 


2  , 


r^^{x,r,t)  dr. 


(19) 


Next,  the  expression  for  is  differentiated  and  the  result  used  to  replace  the  derivative  dcs/dt  in 
Eq.  (19)  in  favor  of  the  derivatives  dT/dt  and  dX/dt,  expressions  for  which  are  given  by  Eqs. 
(13)  and  (14),  respectively.  The  resulting  thermal  conduction  term  vanishes  when  integrated 
due  to  zero  heat  flux  conditions  applied  at  the  boundaries  of  the  energy  localization  sphere  [i.e., 
dT/dr{x,0,t)  =  dT /dr{x,rojt)  —  0].  The  zero  heat  flux  at  implies  symmetry  between  the  uni¬ 
formly  spaced  localization  spheres;  future  work  will  relax  this  condition  to  allow  for  asymmetrical 
thermal  interaction  between  spheres.  Equation  (19)  then  reduces  to 


des 

dt 


3 


2S{x,r,t) 

r  - 

Ps 


dr. 


(20) 
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This  equation  must  be  satisfied  by  our  choice  of  S.  The  thermal  energy  localization  strategy 
specifies  that  bulk  dissipated  energy  is  deposited  within  localization  centers  defined  by  r  <  rc, 
where  the  evolution  of  rc  is  given  by  Eq.  (12).  Therefore,  we  take 


S{x,r,t) 


rcix,t)\ 


3  fjp 


for  0  <  r  <  Tcix^t) 


0 


for  rc{x^t)  <r  <  To 


(21) 


It  is  seen  that  for  S  to  remain  bounded  as  rc  ^  0  requires  that  dcs/dt  rapidly  vanish  in  this  limit. 
While  it  is  not  formally  proven  that  this  condition  exists,  we  note  that  the  singularity  posed  no 
numerical  difficulties  in  this  study. 

In  summary,  the  energy  localization  model  is  given  by  Eqs.  (7),  (12),  (13),  and  (14),  where  S 
is  given  by  Eq.  (21),  Uc  =  and  Vo  -  i?(^7)3.  This  meso-scale  model  is  coupled  to  the  bulk 
model  through  the  thermal  energy  equation,  given  by  Eq.  (5),  and  the  solid  volume  fraction.  The 
composite  model  is  mathematically  closed  and  can  be  solved,  in  principle,  when  suitable  initial  and 
boundary  conditions  are  supplied. 


4  Results 


As  an  illustration  of  the  model,  we  perform  an  analysis  of  localized  heating  and  ignition  of  granular 
HMX  by  steady  compaction  waves  resulting  from  constant  speed  piston  impact.  The  focus  of  this 
analysis  is  on  low  speed  impact  associated  with  weak  initiation  of  DDT  in  heavily  confined  tubes 
of  granular  HMX  [4,  16].  The  analysis  is  similar  to  that  given  by  Gonthier  and  Son  [10]. 


The  model  equations  are  solved  in  a  compaction  wave-attached  frame  defined  by  the  transfor¬ 
mation  ^  =  X  -  Dt  and  v  =  u  -  D,  where  D  is  the  steady  speed  of  a  right  propagating  wave. 
Because  the  model  is  frame  invariant,  the  governing  equations  in  the  wave  frame  have  the  same 
form  as  those  given  in  Section  3  with  x  replaced  by  ^  and  u  replaced  by  v.  The  analysis  is  re¬ 
stricted  to  compaction  waves  propagating  with  speeds  much  less  than  the  ambient  solid  acoustic 
speed  [D  <C  3000  m/s)  due  to  the  incompressibility  eissumption.  We  also  assume  that  PcA/y^c  ^ 
where  Pc  is  a  characteristic  compaction  stress;  thus,  Eq.  (3)  is  replaced  by  the  equilibrium  condition 
^  =  /(^).  We  take  ~  100  kg/(s  m).  For  brevity,  the  constitutive  relations  for  both  /(^)  and 
/3(^,  ^)  are  not  given  here,  but  can  be  found  in  Ref.  [9].  With  these  assumptions,  the  steady  prob¬ 
lem  can  be  reduced  to  the  following  initial-boundary  value  problem  (details  of  the  mathematical 
reductions  are  omitted): 


_  g  d 
vr‘^  dr 


where  a  =  kl{psCy)  and  = 
relation  ^  =  /,  to  obtain  Eq. 


= 

-/?), 

(22) 

vile 

+ 

^  (l  r 

)exp  ^ 

'  T„\  5 

(23) 

dr  j 

<  T  /  ps^'^ 

dX  _ 

(l  -  a)  exp  ( 

-1) 

> 

(24) 

V 

T  J 

drc 

di 

- 

[Ps-p)+r> 

A'Kv'lncPY 
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df  /d(p.  We  have  substituted  Eq.  (5)  into  Eq.  (12),  and  used  the 
(25).  The  bulk  mechanical  response  of  the  granular  material  is 
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given  solely  by  the  solution  of  Eq.  (22)  as  the  variables  v,  P,  jS,  and  ^  depend  only  on  (/>.  Bulk 
temperature  and  reaction  progress  are  related  to  the  meso-scale  response  by  Eqs^  (17)  and  (18), 
respectively.  Initial  conditions  for  the  unstressed  granular  solid  are  <^(0)  =  0.655,  r(0,r)  =  300  K, 
A(0,r)  =  0,  and  rc(0)  =  \.%-nR*Yl{2E*).  Boundary  conditions  for  the  meso-scale  model  are 
df/dr{$,  0)  =  dffdr{^,  ro)  =  0.  The  system  of  partial  differential  equations  (PDEs)  is  numerically 
solved  by  a  Method  of  Lines  technique  using  an  implicit,  stiff  ordinary  differential  equation  (ODE) 
solver  contained  in  the  package  LSODE.  To  this  end,  the  radial  coordinate  of  the  localization  sphere 
is  discretized  into  100  evenly  spaced  nodes  and  the  radial  derivative  in  Eq.  (23)  is  approximated 
by  a  second-order  accurate  centered  finite-difference  relation.  The  numerical  method  was  validated 
against  analytical  solutions  for  compaction  end  states  and  simple  heat  conduction  problems. 

The  bulk  compaction  problem  has  already  been  discussed  in  Subsection  3.1,  and  is  discussed  in 
detail  in  Ref.  [9];  thus,  the  predicted  meso-scale  response  is  the  focus  of  this  section.  We  first  give 
predictions  for  the  localized  heating  of  inert  HMX  (i.e.,  for  A  =  0)  for  Up  =  85.6  m/s,  and  then 
consider  reactive  HMX  for  this  same  value  of  Up.  This  piston  speed  is  chosen  because  it  is  close 
to  that  needed  for  the  predicted  onset  of  sustained  combustion.  The  corresponding  compaction 
wave  speed  is  D  =  297.5  m/s.  The  values  R  =  25  fxm  and  7  =  12  are  used  as  representative  of 
lightly  compacted  granular  HMX.  Model  sensitivity  to  variations  in  these  and  other  parameters 
is  demonstrated.  The  evolution  of  </>,  u,  and  P  within  the  compaction  zone  for  both  of  these 
simulations  are  given  in  Fig.  2. 

All  numerical  simulations  performed  in  this  study  were  solved  on  a  Dell  OptiPlex  GX400  work¬ 
station  having  a  single  INTEL  Pentium  IV,  1.3  GHz  processor.  The  average  computational  run 
time  for  a  single  simulation  was  approximately  5  seconds. 

4.1  Inert  Compaction 

Figure  4  gives  predictions  of  the  energy  localization  model  for  heating  within  the  compaction 
zone  of  inert  HMX.  Grain  number  density  smoothly  increases  from  its  initial  value  of  n  =  1.11  x 
10^  grains/cm^  to  an  approximate  equilibrium  value  of  1.41  x  10^  grains/cm^  through  the  wave 
as  the  material  is  compacted.  The  equilibrium  value  for  the  number  of  energy  localization  sites 
is  Tic  =  =  8.46  X  10^  sites/cm®.  Compaction  induces  mechanical  energy  dissipation  within 

material  contained  by  localization  centers  that  individually  surround  grain  contact  points.  The 
spherical  volume  of  material  contained  within  a  localization  center,  characterized  by  its  radius  r c, 
increases  from  its  initial  value  of  Tc  —  0.93  jjLTCi  to  a  final  value  of  2.55  //m.  As  mentioned,  this 
volume  increase  is  interpreted  as  the  uncontained  spreading  of  plastically  deformed  material  due  to 
intergranular  contact.  The  outer  radius  of  the  localization  sphere  remains  constant  at  to  =  13.76  /xm 
because  the  solid  grains  are  assumed  incompressible.  A  pulse-like  volumetric  heat  deposition  rate 
is  predicted  for  the  locahzation  center,  reaching  the  maximum  value  S  —  0.52  GW/cm^  within 
the  compaction  zone  where  <j>  =  0.73  at  ^  =  -3.71  mm.  Though  this  rate  is  high,  it  is  well  below 
the  characteristic  heat  generation  rate  of  approximately  Sd  =  psQ^dlLd  =  3.1  x  10^  GW/cm^  (for 
Dd  —  7000  m/s  and  Ld  =  2.55  /am)  that  would  be  associated  with  detonation  of  the  solid  HMX 
contained  within  the  localization  center.  This  heating  gives  rise  to  the  temperature  field  of  Fig. 
4(e).  The  plot  shows  radial  temperature  profiles  within  a  localization  sphere  though  the  compaction 
zone.  A  maximum  temperature  of  735.5  K  is  predicted  at  the  center  of  the  localization  sphere  at 
^  =  —4.69  mm;  thus,  the  peak  temperature  within  the  compaction  zone  trails  the  peak  energy 
deposition  rate.  The  temperature  rapidly  decreases  as  the  applied  energy  deposition  rate  vanishes. 
Though  not  shown  here,  the  temperature  field  within  the  localization  sphere  will  approach  the 
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uniform  bulk  value  of  T  =  302.4  K  as  ^  ^  -oo.  This  feature  is  a  consequence  of  having  energetically 
consistent  localization  and  bulk  models. 

Inspection  of  Figs.  4(c)  and  (e)  indicates  that  much  of  the  applied  thermal  energy  remains  in  the 
vicinity  of  the  localization  center  through  the  compaction  wave  as  there  is  little  time  for  conduction 
into  the  cooler  grain  interior.  This  observation  does  not  imply  that  thermal  conduction  is  unim¬ 
portant  in  the  heating  process.  An  estimate  for  the  change  in  temperature  of  material  contained 
in  the  localization  center  in  the  absence  of  thermal  conduction  is  given  by  AT  =  SavgT/ipsCv), 
where  Savg  and  T  are  the  average  applied  energy  deposition  rate  and  grain  residence  time  within 
the  compaction  zone,  respectively.  Taking  Savg  —  0-2  GW/cm^  and  T  =  19  ps,  which  are  charac¬ 
teristic  of  the  predictions  shown  here,  then  AT  =  1331  K.  Indeed,  the  same  simulation  performed 
with  a  =  0  m^/s  predicts  a  temperature  within  the  localization  center  that  is  in  excess  of  1400  K. 
We  further  determine  the  mass  fraction  of  heated  material,  X,  having  temperature  within  intervals 
of  approximately  10  K  over  the  range  300  <  T  <  750  K;  the  result  is  shown  in  Fig.  4(f).  Only 
a  small  fraction  of  the  total  mass  is  heated  to  elevated  temperatures  within  the  compaction  zone, 
and  this  mass  is  rapidly  cooled  by  thermal  conduction  as  the  energy  deposition  rate  diminishes.  An 
important  question  is  whether  the  heated  mass  has  sufficient  thermal  inertia  to  trigger  the  onset 
of  sustained  combustion.  This  question  is  answered  in  the  following  subsection. 

We  lastly  remark  that  phase  change,  not  accounted  for  by  the  model,  will  affect  predicted 
temperatures.  A  relevant  discussion  is  given  by  Menikoff  and  Kober  [17].  While  phase  change 
would  reduce  predicted  temperatures  due  to  the  latent  heat  of  fusion,  thermal  energy  localization 
in  real  materials  will  likely  be  more  pronounced  than  predicted  here  due  to  stress  chain  formation. 
As  such,  these  predictions  are  considered  a  reasonable  average  localization  response.  The  inclusion 
of  phase  change  will  be  addressed  at  a  later  time. 

4.2  Reactive  Compaction 

Figure  5  gives  predictions  of  the  energy  localization  model  for  heating  and  chemical  reaction  within 
the  compaction  zone  of  HMX.  Profiles  for  n,  ric,  Vc,  To,  and  S  are  the  same  as  those  given  in  Fig. 
4.  The  predicted  variation  through  the  compaction  zone  of  the  radial  temperature  profile  within 
a  localization  sphere  is  shown  in  figure  (a);  the  corresponding  variation  in  local  reaction  progress 
is  shown  in  (b).  The  temperature  near  the  center  of  the  sphere  is  sufficient  to  induce  chemical 
reaction  that  locally  goes  to  completion  (A  =  1).  The  reaction  accelerates  almost  to  the  point  of 
thermal  explosion  as  indicated  by  the  presence  of  a  sharp  temperature  spike.  However,  thermal 
energy  conduction  from  the  reaction  site  into  the  grain  interior  freezes  the  reaction  causing  a  rapid 
decrease  in  temperature.  A  peak  temperature  in  excess  of  1200  K  is  predicted  prior  to  extinction 
of  chemical  reaction.  The  predicted  evolution  of  hot-spot  mass  fraction  is  shown  in  figure  (c).  As 
before,  only  a  small  fi-action  of  the  total  mass  is  heated  to  elevated  temperatures,  and  its  thermal 
inertia  is  insufficient  to  overcome  thermal  conduction  losses  and  achieve  sustained  combustion. 

The  predicted  variation  in  bulk  temperature  and  reaction  progress  through  the  compaction 
zone  are  shown  in  Fig.  6.  Because  bulk  temperature  depends  on  the  temperature  field  within  a 
localization  sphere  by  Eq.  (17),  the  profile  is  not  identical  to  that  that  given  in  Fig.  2(d)  for  the 
inert  case.  However,  since  chemical  reaction  is  mostly  confined  to  a  very  small  volume  of  material 
near  the  sphere  center,  the  equilibrium  value  of  bulk  temperature  following  compaction  is  also 
approximately  T  =  302.4  K. 

For  this  simulation,  Up  =  85.6  m/s.  Numerical  experiments  indicate  that  this  piston  speed  is 
close  to  the  thermal  explosion  threshold;  the  value  Up  =  85.7  m/s  results  in  a  thermal  explosion 
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Figure  4:  Predicted  localized  heating  through  the  compaction  zone  for  inert  HMX  {up  =  85.6  m/s): 
(a)  grain  number  density,  (b)  localization  center  number  density,  (c)  localization  center  radius,  (d) 
volumetric  heating  rate,  (e)  local  temperature,  and  (f)  hot-spot  mass  fraction. 
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marking  the  onset  of  sustained  combustion  of  the  material.  This  prediction  agrees  well  with  the 
results  of  confined  DDT  tube  tests  involving  piston  impact  of  granular  HMX  [16,  22].  These 
tests  indicate  that  the  explosion  is  first  triggered  for  piston  speeds  approximately  in  the  range 
80  <  Up  <  90  m/s.  Combustion  is  often  observed  for  piston  speeds  slightly  below  this  threshold, 
but  transition  to  detonation  does  not  occur. 


T[K]  ^ 


(a)  (b) 


X 


Figure  5:  Predicted  heating  and  reaction  progress  through  the  compaction  zone  for  HMX  {up  — 
85.6  m/s):  (a)  local  temperature,  (b)  local  reaction  progress,  and  (c)  hot-spot  mass  fraction. 


4.3  Parametric  Sensitivity 

Numerical  experiments  were  performed  to  determine  the  model  sensitivity  to  variations  in  key 
energy  localization  parameters.  To  this  end,  parameter  values  used  for  the  reactive  compaction 
simulation  were  chosen  as  a  baseline,  and  the  values  of  7,  PcIY,  and  Py/Y  were  independently 
varied  over  physically  significant  ranges.  The  piston  speed  needed  for  sustained  combustion  (ther- 
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(a)  (b) 

Figure  6:  Predicted  bulk  (a)  temperature  and  (b)  reaction  progress  through  the  compaction  zone 
forHMX  =  85.6  m/s). 

mal  explosion)  was  determined  for  each  case.  Here,  the  value  of  7  determines  the  hot-spot  density 
of  the  meso-structure.  The  ratio  PdY  determines  the  initial  radius  of  the  localization  center,  rc(0), 
as  seen  from  Eq.  (10),  and  the  ratio  Py/Y  determines  the  growth  rate  of  given  by  Eq.  (12). 
Thus,  both  of  these  ratios  control  the  volume  over  which  bulk  dissipated  energy  is  localized.  Be¬ 
cause  grain  size  is  an  important,  easily  adjustable  system  parameter,  its  influence  on  the  predicted 
onset  of  sustained  combustion  was  also  determined. 

Results  of  the  sensitivity  analysis  are  summarized  in  Fig.  7.  The  value  of  7  was  varied  over 
the  range  1  <  7  <  14;  it  is  reasonable  to  expect  that,  on  average,  the  number  of  contact  points 
per  grain  lies  within  this  range.  The  predicted  threshold  for  sustained  combustion  increases  with 
7  as  bulk  dissipated  energy  is  distributed  over  a  larger  material  volume.  Values  for  Pc/Y  and 
Py/Y  were  varied  over  the  ranges  0.1  <  Pc/Y  <  3.0  and  1.6  <  Py/Y  <  6.0  which  are  typical 
of  contact  problems  [13].  The  sustained  combustion  threshold  is  predicted  to  increase  with  Pc/Y 
and  decrease  with  Py/Y]  these  results  are  a  consequence  of  bulk  dissipated  energy  being  localized 
over  larger  and  smaller  material  volumes,  respectively.  Importantly,  all  of  these  predictions  are 
near  the  experimentally  observed  explosion  threshold  for  granular  HMX.  This  result  suggests  that 
the  localization  strategy  used  in  this  study  may  be  a  viable  approach  to  develop  improved  ignition 
models  for  heterogeneous  reactive  solids. 

Lastly,  the  influence  of  grain  size  on  the  predicted  combustion  threshold  is  shown  in  Fig.  7(d). 
Here,  it  is  seen  that  the  piston  speed  for  sustained  combustion  rapidly  increases  as  grain  size 
decreases.  A  value  of  Up  =  384.2  m/s  is  needed  to  trigger  sustained  combustion  for  a  material 
composed  of  R  =  2.5  ^m  size  grains.  This  result  agrees  with  experimental  observations  which 
show  the  same  trend.  Further,  the  necessary  piston  speed  decreases  only  slightly  for  grains  having 
approximately  i?  >  30  fj,m. 

5  Conclusions 

This  report  summarized  a  model  and  analysis  of  mechanically  induced  hot-spot  formation  and 
ignition  of  granular  HMX  that  is  associated  with  weak  initiation  of  DDT.  To  this  end,  a  bulk 
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Figure  7:  Sensitivity  analysis  results:  (a)  contacts/grain,  7;  (b)  initial  localization  center  radius, 
characterized  by  Pc/^;  (c)  localization  center  radius  growth  rate,  characterized  by  Pyiy ;  and  (d) 
grain  size,  R. 


compaction  model  was  coupled  with  a  meso-scale  model  that  tracked  the  evolution  of  hot-spots 
within  the  material  meso-structure.  Importantly,  the  interplay  between  localized  energy  deposition, 
thermal  conduction,  and  chemical  reaction  at  the  grain  scale  was  resolved.  The  deposition  of  bulk 
dissipated  mechanical  energy  in  the  vicinity  of  intergranular  contact  points  was  guided  by  basic 
principles  of  contact  mechanics.  The  bulk  and  meso-scale  models  were  coupled  in  an  energetically 
consistent  manner  which  preserved  the  integrity  of  the  experimentally  correlated  bulk  predictions. 

Localized  heating  predictions  indicate  that  temperatures  sufficient  to  induce  chemical  reaction 
can  result  from  weak  piston  impact.  Temperatures  near  700  K  were  predicted  for  piston  speeds 
near  86  m/s,  the  value  needed  for  the  onset  of  sustained  combustion.  This  result  agrees  well  with 
experimentally  determined  explosion  thresholds  for  granular  HMX.  Also,  thermal  conduction  was 
shown  to  have  a  significant  influence  on  the  ignition  process;  in  the  absence  of  thermal  conduction, 
excessively  high  temperatures  were  predicted  near  intergranular  contact  points.  The  inclusion 
of  thermal  conduction  is,  thus,  believed  essential  for  accurately  describing  the  weak  initiation  of 
DDT.  Results  of  this  study  suggest  that  mechanically  induced  ignition  can  be  accurately  predicted 
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if  important  meso-scale  features  responsible  for  hot-spot  formation  are  resolved. 

Predictions  showing  the  model  response  to  variations  in  key  energy  localization  parameters 
indicate  that  it  is  reasonably  insensitive  to  the  localization  strategy.  Because  the  model  already  gives 
results  that  agree  well  with  experiments,  it  is  possible  that  very  good  predictability  can  be  obtained 
by  including  a  more  detailed  description  of  the  localization  process  including  compressibility,  stress 
phaiu  formation,  phase  change,  and  grain  fracture.  Such  a  description  would  greatly  increase  model 
complexity  and  thus  necessitate  the  use  of  modern  computational  tools  for  its  accurate,  efficient, 
and  timely  solution.  Information  provided  by  meso-scale  experiments  and  numerical  simulations 
can  be  used  to  guide  the  development  of  improved  energy  localization  strategies. 
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